Stellar Atmospheres Near an AGN: 
The Importance of Radiation Pressure from Trapped Lyman-a Photons 
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ABSTRACT 



We derive an analytic expression for the intensity of resonance-line radiation 
Q\ ' "trapped" in a semi-infinite medium. Given a source function and destruction 

o ' 

(T) ■ calculated by numerically integrating over analytic functions. We apply this formalism 



to a plane-parallel model stellar atmosphere to calculate the radiation pressure due 
to Lyman-a photons produced following absorption of UV and X-rays from an AGN. 
For low surface gravity stars near the AGN (g ~ 10 cm s~ 2 , r ~ 0.25 pc), we find that 
the pressure due to Lyman-a photons becomes an appreciable fraction of that required 
for hydrostatic support. If the broad emission line emitting gas in AGNs and QSOs 
' consists of stellar outflows, it may be driven, in part, by Lyman-a pressure. 



Subject headings: galaxies: active — line: formation — radiative transfer — stars: 
atmospheres 



1. Introduction 

Interest in the effect of hard UV and X-ray radiation on stellar atmospheres and stellar 



evolution began in earnest with work on X-ray binaries (e.g., [Davidson fc Ostriker 1973| , |Basko fc 



Sunyaev 1973| , and |Arons 1973Q , but the subject has since been revisited in the context of active 



galactic nuclei (AGN). 

Consider a star of effective temperature T e g- and an AGN of luminosity L = 10 46 L^q erg s . 
The incident power per area from the AGN is equal to o~T^ s at the "heating" distance 
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heating of the stellar surface by the AGN is important at distances d <^ dh- Fabian (1979) was 
apparently the first to suggest that heating by AGN radiation might substantially affect the 
envelopes of nearby stars, resulting in enhanced mass loss. Edwards (1980) argued that when a 
star approached within d <^ dh, the irradiated photosphere would develop supersonic horizontal 
winds carrying heat from the illuminated hemisphere to the shadowed side. Edwards conjectured 
that enhanced mass loss could then occur, resulting in a "cometary star" as the stellar outflow 
is accelerated radially away from the AGN. Matthews (1983) noted the importance of direct 
radiation pressure from the AGN, arguing that this could ablate matter tangentially from the 
stellar photospheres of giant or supergiant stars. 

Voit & Shull (1988) discussed the effects of X-rays on the upper atmosphere of a star near 
the AGN. They argued that X-ray heating of a ~ IMq, i?* ~ IOORq star (g ~ 3 cm s -2 ) would 
result in a hot, ionized wind, with temperature at the critical point T ~ 3 x 10 5 K; radiation 
pressure in C IV, N V, and O VI resonance lines would contribute to the acceleration of the 
wind. They concluded that stars which are already red supergiants could develop winds with 
M « lO~ 7 L^ 9 dj" 8 2 i? 2 iloo M yr- 1 , where the distance from the AGN is d = 10 18 di 8 cm, and 
R* t ioo = R*/W0Rq. These mass loss rates exceed normal mass loss rates for red supergiants 
only within a distance of d\s <C, 0.3L^q 5 R^^qq cm. They estimated the mass loss due to ablation 
by radiation pressure dMatthews 1983|) for a ~ 1M & star to be M « lO^L^d^ 2 ^^ provided 



y 2 

d\s < 0.6L 46 -R*,ioo- 

There has also been some exploration of how AGN radiation might affect the evolution of 
stars ( [Matthews 1983 ; Verbunt, Fabian, &; Rees 1984; Tout et al. 1989| ). Tout et al. discussed 



the evolution of stars immersed in a blackbody background with temperature T = 10 3,75 K, 
concluding that the main sequence evolution is largely unaffected, but that the star will expand 
to a radius ~ 10 times larger than usual during the "red giant" phase of evolution. Note, 
however, that the energy density of AGN radiation equals a 10 4 K blackbody at a distance of only 
2 x 10 16 L4g 5 cm from the AGN. Furthermore, it is not clear how their conclusions would have to 
be modified for the anisotropic and nonthermal radiation field of the AGN. 

These discussions of the effects of AGN radiation on stellar atmospheres and evolution have 
led some to consider a "bloated stars scenario" (BSS) to explain the origin of the AGN broad 
line region (BLR). The BSS proposes that the BLR consists of the stellar winds and/or expanded 
envelopes of stars irradiated with AGN radiation (e.g., Penston 1988 , Kazanas 1989j ), The 



attractiveness of this scenario stems from its efficient use of known resources: stars are believed 
to be present near the "central engines" of AGN, their gravity provides a possible "containment" 
mechanism for BLR clouds, there is ample mass with which to replenish the clouds, and the 
mass shed provides material to fuel the AGN. So far, however, no specific model has successfully 
explained how all these mechanisms might function, though pieces of the puzzle have been explored 
(Alexander & Netzer 1994). High signal/noise measurements of line profiles indicate that the 



number of discrete clouds or bloated stars must be large, perhaps exceeding ~ 3 x 10 6 for Mrk 335 
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(Arav et al. 1997) and ~ 3 x 10 7 for NGC 4151 (Arav et al. 1998). Such large numbers appear to 
challenge the BSS, but the difficulties may not be insurmountable. 

Recent observations and photoionization models of Baldwin et al. (1996) and Ferland et al. 
(1996) offer some concrete evidence for the BSS. They find that the BLR seems to contain several 
distinct components. One (component "A") has sharp (FWHM ~ 1000 km s -1 ), symmetric line 
profiles centered on zero velocity, while another (component "B") appears to be outward-flowing, 
peaked at zero velocity but with a long blue tail (down to —11,000 kms -1 ). They interpret "A" 
as the expanded envelopes of stars, and "B" as their radiatively accelerated, outflowing winds. 

The aim of the present paper is to call attention to the fact that irradiation of a star by hard 
X-rays from the AGN will result in the generation of Lyman-a (Lya) photons within the stellar 
atmosphere, and the pressure of these trapped resonance-line photons may contribute to mass loss. 
The possible importance of Lya photons was suggested previously by Puetter [unpublished, cited 
in Penston (1988) and Voit & Shull (1988)]. Voit & Shull rejected Puetter's suggestion, arguing 
that the Lya pressure could be estimated to be 



where W\ is the equivalent width of the Lya absorption profile of the gas between the point of 
interest and the surface, and T is the gas temperature. Voit &: Shull then used W\/\ < 1 to obtain 
an upper limit Phja ^ aT 4 /3 which, for stars of interest, is insufficient to "bloat" the atmosphere. 



the gas temperature T — an assumption which need not be valid when Lya photons are being 
generated within the gas by nonthermal processes, such as H(ls— >2p) excitation by photoelectrons 
and secondary electrons in a non-Maxwellian "tail" to the electron energy distribution function. 
It is therefore essential to estimate the rate of production of Lya photons, and to examine their 
diffusion in physical space and frequency space as well as the possibility of photon "destruction" 
by, for instance, collisional de-excitation of electronically excited H atoms. 

In §|2[ we give a brief treatment of resonance- line transfer, extending the detailed calculations 
of Neufeld (1990). In particular, we derive an expression for the case of a semi-infinite, absorbing 
slab of material. Since resonance-line photons created in a stellar atmosphere will almost certainly 
not be able to scatter through the entire star without being absorbed, this limit is an excellent 
approximation. In §^, we apply our calculational method to a simple, plane parallel, irradiated 
model atmosphere. In §|4], we present the results. In §|5] we summarize and discuss implications 
and further directions for work. 




(2) 




2. 



Resonance Radiation Transfer 



Most current photoionization treatments, including CLOUDY (Ferland 1996), use the escape 
probability approximation to address resonance-line radiation transfer. For example, Elitzur and 
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Ferland (1986) use the escape probability approximation (EPA) to derive the frequency-integrated 
radiation intensity J u l, EPA for a resonance line at frequency v u i, corresponding to a transition 
from level u to level I, to be 

Jul, epa = B(v u i, T exc )Az/ M /, (3) 

where B(iy u i,T exc ) is the Planck function at the line excitation temperature (which is calculated 
from the level populations) and Av u i is the line width. The radiation pressure P u i, epa is given, as 
usual, by 

4ir 

Pul, EPA = -^-Jul, EPA- (4) 
6C 

In this purely local treatment, the escape probability approximation assumes that a resonance 
photon at a point r has a probability p u i per scattering of "escaping" to a defined boundary 
ro, where it is assumed photons free-stream to infinity. Thus, the differential radiation pressure 
P^ lfl EPA (r) at any point r is given by 

*5epaW = ^,EPA(r)-P^ E PA(ro) (5) 

= Pul,EPA(r) x {I -Pul) (6) 
4-7T 

= ~^B(v uh T cxc )Av ul (l -pui). (7) 

Confusingly, it is this differential radiation pressure which Elitzur and Ferland (1986) and 
CLOUDY (Ferland 1996) refer to as the "line radiation pressure." We will always use the term 
"pressure" and the symbol Pi for the physical pressure (momentum flux) at point r due to particles 
i, such as gas particles or resonance line photons. We will use the term "differential pressure" and 
the symbol Pf lS = Pi(r) — P«(ro) for the pressure difference due to particles i between the point 
of interest r and the boundary point r . 

Elitzur and Ferland (1986) note that the line width Av u i is the "crucial parameter" in this 
calculation of the radiation pressure. In particular, it is not well understood how to calculate Av u i 
at extremely large optical depths when the escape probability is comparable to the "destruction 
probability" per scattering, for example due to collisional de-excitation. 

However, the case of very large optical depth is precisely that which Neufeld (1990) considers 
using a completely different approach. He presents in §§2-3 of his paper an exact and analytic 
calculation derived from Fourier transforming a diffusion equation. He solves this equation for 
a unit delta-function photon source in a uniform, optically thick (but finite) slab, with and 
without a finite destruction probability. In particular, at optical depth r with a source at t s , 
the angle-averaged intensity J(t,t s ,x) and frequency-integrated intensity J(t,t s ) are given in 
terms of infinite series [e.g., his equation (2.7)]. Unfortunately, these series converge rather slowly. 
However, we show below that if we take the limit of a semi- infinite slab, then the angle- averaged 
intensity J(t,t s ,x) can be expressed in terms of analytic functions which can be integrated 
numerically with relative ease. Below, we first present the notational conventions we use, and then 
proceed with our derivation of the intensities J and J in the case of a semi-infinite slab. 
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Neufeld (1990) defines a as the ratio of the natural line width to the Doppler width and r as 
the mean optical depth:[] 

7 A _ vre 2 //«A , 
r = a a N , a a = -= (8) 



where 7 is the decay rate, Ao is the line-center wavelength, = \JkT '/ran is the one-dimensional 
velocity dispersion, mn is the mass of atomic hydrogen, fi u is the oscillator strength, m e is the 
electron mass, Nq is the column density of absorbers, and the other symbols have their usual 

meanings. For Lya, this corresponds to v\d = 9.082 T 4 kms -1 , a = 4.717 x 10~ 4 T 4 , and 
' _ 1/ 

r = 1.045 x 10~ 13 T 4 Nq cm 2 , where T4 is the temperature in units of 10 4 K. The dimensionless 
frequency variables x and cr(x) which Neufeld uses are given by: 

~- Xo{u - uo) *( X) =Jirj*, ( 9 ) 



V2vi d ' V 3 Jo 4>{x') 

where 4>(x) is the Voigt profile normalized such that / <ft(x)dx = 1. In this convention, (j)(x) takes 
on the limits 

0(0) = — (1 + 0{a)) ; 0(a;)«-^,x»l. (10) 
'it ttx z 



The intensity J(t,t s ,x) is given per unit frequency in x, so that the frequency-integrated intensity 
is given by 

J(t,t s ) = J J(r,r s ,x)dx. (11) 

Since Neufcld's solutions for the intensities J{t^t s ^x ) and J(t,t s ) are derived for a 
delta-function source, they act as Green's functions for resonant photon sources. Consider a slab 
of total thickness To- For a Lya emissivity per unit volume 47rj, we define the injected power per 
unit area per unit optical depth as 

«r.) = *gp>. (12) 

The excitation temperature T exc at optical depth r is found by integrating J near line center over 
all source locations r s , weighted by £, to obtain the total intensity at line center, and comparing 
with a blackbody intensity: 

B(v ,T exc )-^ = J o Z(r s )J(T,T s ,0)dT s , (13) 

where -B(fo,T exc ) is the usual Planck function at the excitation temperature T exc , evaluated at 
line center. Since the majority of scatterings occur near line center, the local level population for 
Lya should be given by 

^-^expf-^V (14) 



For consistency, we retain the definitions of r and a a used by Neufeld (1990), Harrington (1973), and other 
workers: a a — \/na , where a is the absorption cross section at line-center. 
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where the n's are the level populations, g's are the statistical weights, and X21 = 10.2 eV. The 
radiation pressure P tiu i(t) = 47r l 7/(3c) is given similarly by 



Air f T ° 

J Prad(r) = — / £(T s )J(T,T s )dT s . 
3c Jo 



(15) 



Consider photons created at line center with a constant probability per scattering eo <C 1 that 
the photon will be "destroyed" by, for instance, collisional de-excitation. Using Neufeld's equations 
(2.7) and (3.11)^(3.13), we find that for a slab of total thickness tq and source at < r s < tq, the 
angle-averaged intensity at < r < tq can be written as 



J(t,t s ,x) 



8vr 2 I 



F[t — T s ,a(x)\ — F 



t + t s + 



+ O[r 1 ),a(x 



F[v,a{x)} = — ]T COS ^ n ^,„ exp[-/i ra |cr(x)|], 



Mr! 



To n=1 [l n + eoA/6/2 

— i + ofv 1 

r L v 



(16) 
(17) 
(18) 



Note that our tq is twice Neufeld's tq. 



In a stellar atmosphere, we essentially have a semi-infinite slab, where to — > oo. We can 
therefore neglect the terms of order r _1 , and replace the sum in equation ( |17| ) with an integral. 
With the substitutions 

Ta = 75, (19) 



v/6' 



T d [i n = nirT d /To, —— = VTT d /r , 
an 



F then becomes 



F(v,a) 



cos(uv/Td) 



z = 



1 + u 

£(z) = Re[e^i(z)] 
\<j(x)\ + i\v\ 



exp[— u\a(x) | /Td]du 



Td 



(20) 

(21) 
(22) 
(23) 



where Ei(z) = J 2 °° w l e w dw is the exponential integral with complex argument. The quantity Td 
is the effective optical depth due to finite destruction probability. 



Therefore, the angle-averaged intensity in the limit that tq — > oo is given by 



J(t,t s ,x) = ^{8 



\a(x)\ + i\r - t s 

Td 



\a(x)\ + i\r + t s + 



3<t>(x) i 



Td 



(24) 



At line center, where x = a = and the argument z is purely imaginary, £{z) can be expressed in 
terms of cosine and sine integrals Ci(y) = — f°° t~ l cos t dt and si(y) = — t~ l sin t dt: 

£(iy) = ~ [cos y Ci(y) + sin y si(y)) . (25) 
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Using equations (|24|)— (p5|) in equation (13) gives the excitation temperature, and hence the level 
population at each point in the slab using equation (|l~4|). 

To find the frequency-integrated intensity J , we must integrate J over the frequency x, 
equation flll]). The exact expression for the frequency- integrated intensity in the limit that 
ro — ► oo is thus: 

r°° \cr{x) + iv 

10 [ T d 



= 2 x 



Gr d (r + T s 

dx, 



(26) 
(27) 



2 x 



£ 



<7{X 



)+<( 



dx. 



(28) 



The functions J- and in general need to be calculated numerically. Obtaining the numerical 
integral is simplified by taking into account the following limits of the integrand: 



£ 



<t(x) + iv 



constant if x ^ (av) 1 3 and x <J (ar d ) 1 3 ; 



oc < ~ X £ 



~ X 



if (ar d ) / 3 ^ x £ (aw) / 3 and T d < w; 
if x ^ (ar^) / 3 and x ^ (aw) / 3 . 



(29) 



If the Lya optical depths r and r s , as well as the effective optical depth to destruction r^, 
are all large, so that (av) / 3 3> 1 and (ar d ) / 3 3> 1, then the scattering is dominated by the 
damping wings. We may then approximate the Voigt profile <f>{x) by using equation (10) and <r(x) 
by cr(x) rj y/2/3 7rx 3 /(3a). This approximation leads to an analytic expression for J 7 and £/: 



Or d (v) 



( 1/3) r ( 2 / 3 ) f -^-) 3 Re [exp(^/r d )r ( Vs, iv/r d )] 



\V2 



7T 



ft. 



V + 



3<A(0) 



(30) 
(31) 



where the incomplete gamma function with an imaginary arg ument is F(b, iy) = i b f™ t^e^dt. 
This case is not applicable to the AGN-illuminated stellar atmosphere model we discuss below, 
but we include it for completeness. 



3. A Stellar Atmosphere Near an AGN 



Our procedure for evaluating the importance of Lya pressure in a stellar atmosphere near an 
AGN is as follows: 



A. Use a slightly- modified version of CLOUDY (Ferland 1996) to calculate the thermal and 
ionization structure of a plane-parallel, gravitationally confined gas illuminated by an AGN. 
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B. Calculate the Lya production and destruction rates from A. 

C. Use B and the analytic model from §2 to calculate the Lya radiation pressure Ph ya ( r )- 

D. Compare the differential radiation pressure P^ya = PLyaM — Phya( r o) from C with Ph y d ; 
the pressure required for hydrostatic equilibrium. 

If ^LyL ~ then the atmosphere may be unstable, especially since other resonance lines will 

likely contribute to the total radiation pressure (Elitzur & Ferland 1986). We describe each of 
these steps in greater detail below. 

3.1. Use of CLOUDY 

We have modified the CLOUDY program (version 90.02) of Ferland (1996) to model a 
plane-parallel atmosphere with a constant gravitational acceleration, aligned normal to the 
incident AGN flux. We ignore the radiation from the star, considering only the effect of the AGN 
illumination on a gravitationally confined gas. The modifications to CLOUDY thus only involve 
adding a constant gravitational acceleration g to the equation for the total pressure equilibrium: 

[OcontM + g]p(r)dr + Pgas(r ) + PlinesM- (32) 

Here r is the distance from the center of the star, ro is the starting point of the integration 
(r < ro), -Pgas(^) and Pines (0 are the pressures due to gas and resonance lines respectively (the 
latter calculated by CLOUDY using escape probabilities), and p{r) is the gas mass density. 
We define the hydrostatic pressure P ny d to be the total pressure P to t minus the free-streaming 
radiation pressure at ro, Pii ne s( r o) : 

Phyd(r) = Ptot(r) - Pines (ro). (33) 

This hydrostatic pressure provides the basis of comparison with the differential radiation pressure 
Shies' since it i s the equivalent gas pressure necessary for hydrostatic equilibrium, in the absence 
of any line radiation pressure. The acceleration due to the momentum deposition from the 
attenuation of the AGN continuum flux is given by 

1 [°° s 
Ocont(^) = / <?l (r)fv(r)dv, 

where is the gas mass per H nucleon in units of the mass of atomic hydrogen mu, fv{r) is 
the AGN energy flux per unit frequency, and <r£ ff (r) is the effective cross section per H nucleon, 
the latter two quantities evaluated at frequency v and distance r. Recall that the line radiation 
pressure calculated by CLOUDY Plwj OUDY (?') is actually the differential line radiation pressure 
^L?ya( r )- Our treatment, on the other hand, explicitly calculates both PLya(^) and PLya(^o)) 



(34) 
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the momentum flux in free-streaming resonance photons at ro, and takes their difference to find 



The abundances we use are those used by Ferland et al. 1996 and Baldwin et al. 1996 



log 10 (H: He: Li: Be: B: C: N: O: F: Ne: Na: Mg: Al: Si: P: S: CI: Ar: K: Ca: Sc: Ti: V: Cr: 
Mn: Fe: Co: Ni: Cu: Zn) = (0: -0.907: -7.92: -9.82: -8.35: -3.08: -2.94: -2.23: -6.71: -3.10: 
-4.87: -3.6: -4.71: -3.64: -5.62: -3.97: -5.92: -4.63: -6.06: -4.85: -8.76: -6.91: -7.83: -6.16: 
-6.31: -4.18: -8.5: -5.6: -7.58: -7.19). The incident AGN flux we use is the same one used by 



Baldwin et al. 1996 and Ferland et al. 1996 in their photoionization models of the BLR. Their 



formula is equivalent to 

- lb = M%-"" +0M '-*") (35> 

hv < £ max 

where hv c = 21.5 eV and hv® = 13.6 eV, which corresponds to a UV-X-ray spectral index 
a ox = —1.2. For a Hubble's constant Hq = 65 /i65 kms -1 Mpc , Baldwin et al. (1996) adopt 
A = 1.77 x 10 47 /ig 5 2 erg s _1 and a distance of d = 0.25 pc, both of which we adopt as well. 
We use a high energy cutoff of E max = 100 keV. The luminosity in ionizing radiation is related to 
equation (j35"|) by 

Lion = / L v dv = 1.57 • A. (36) 

Using Lj on as an estimate for the AGN luminosity heating of the star, we obtain using equation 
(0) a "heating distance" of 

1/ 



=02 H^J' ! - PC ' (37> 

Therefore, for a T e ^ = 5000 K star, our calculations are at the heating distance. 

Finally, a few notes about our integration procedure in CLOUDY. We begin our integration 
at a large distance with a total hydrogen density of 10 10 cm -3 [which also determines -P ga s(^o) 
in equation fl32|) 1 and an unattenuated incident flux (p5|); we integrate inward and stop when 
the electron fraction drops to 5%; and we iterate the solution until pressure equilibrium (with 
CLOUDY's implementation of the resonance-line radiation pressure) is reached. 



3.2. Calculation of Lyman-a Production and Destruction Rates 

Given the ionization equilibrium calculated by CLOUDY, we calculate the Lya production 
and destruction rates in the following manner. For Lya production, we assume all electrons 
entering either n = 2 states which do not involve the absorption of a scattered Lya photon result 
in the production of a new Lya photon. We are justified in including electrons entering the 
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2s state, because at the densities we consider, an electron in the 2s state has an overwhelming 
probability of being collisionally transferred to the 2p state, and because once an electron is in the 
2p state, the overwhelming probability is that it will result in the emission of a Lya photon. Thus 
the production processes we consider include (1) radiative and collisional (3-body) recombination 
directly to n = 2, (2) photoexcitation by the AGN continuum, (3) collisional excitation by thermal 
electrons, (4) radiative and collisional de-excitation from n > 2, and (5) collisional excitation by 
photoelectrons and secondary electrons. This gives a Lya production emissivity per unit volume, 
as a function of hydrogen column density iVn, of 

5 

4vri(Ar H ) = X 2in H E r iH 2p, 2s) (38) 
i=i 

where X21 is 10-2 eV, njj is the total hydrogen number density, Tj is the rate for process j per 
hydrogen nucleon. In Figure 1, we show the Lya photon production rate per hydrogen nucleon [i.e., 
the summation in equation (p8|)] in the case of stellar atmospheres with g = 1, 10, and 100 cm s -2 . 

The destruction probability eo is the probability per scattering that either an electron in 2p 
does not end up producing a (scattered) Lya photon, or that a scattered Lya photon is absorbed 
by something other than an electron in the Is state of hydrogen. Thus the processes we consider 
are (1) photo- and collisional ionization from 2p, (2) collisional de-excitation from 2p to Is, (3) 
photo- and collisional excitation from 2p to n > 2, (4) absorption by the H2 Lyman bands from 
v = 2, J = 5, 6, and (5) transition to 2s followed by (5a) two-photon decay to Is, (5b) photo- 
and collisional ionization from 2s, (5c) collisional de-excitation from 2s to Is, (5d) photo- and 
collisional excitation from 2s to n > 2. For molecular hydrogen, we use the rates and population 
fractions from Neufeld (1980, section V, a, ii-iii). We therefore obtain a destruction probability 
per scattering of: 



e = 1.2 x 1Q-^ + x 
n H 



T(2p 2s) 



rx(2 P -) + r 2 (2 P -) + r 3 (2 P -) + 
r 5a (2s -0 + r 5b (2s -►) + r 5c (2s -►) + r 5d (2s -> 



(39) 



T(2s 2p) 

where the first term is due to molecular hydrogen, A21 is the Lya decay rate, Tj(2p — ►) is the 
rate per 2p atom for process i, Ti(2s — ►) is the rate per 2s atom of process i, and T(2p — * 2s) and 
r(2s — > 2p) are the rates per 2p or 2s atom of transition to 2s or 2p, respectively. In Figure 2, we 
show the Lya destruction probability per scattering which we calculate using equation (^) and 
the results of CLOUDY for g = 1, 10, and 100 cms -2 . We note that we find Lya destruction to 
be dominated by photo- and collisional ionization and excitation out of the n = 2 states. 



4. Results of Calculation of Lyman-a Pressure 



With the distribution of Lya sources from equations (|T 
equation ( |i~5| ) with equations p^)-(p8|) and (| 



and (|38[), we numerically integrate 
to determine Phya i n each layer. 
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For the incident flux Aj{^d 2 ) = 1.77 x 10 47 erg s' 1 /[4tt(0.25 pc) 2 ] estimated by 
Baldwin et al. 1996 and [Ferland et al. 1996| (corresponding to a luminosity L46 = 27.8 /ig 5 2 ), the only 



parameter that can be varied is the gravitational acceleration. We consider g = 1, 10, 100 cm s~ 2 . 
The results are presented in Figures 3-5. It appears that for g <^ 10 cms -2 , the Lya pressure 
rises to ~ 10% of the hydrostatic pressure. In this region, the total hydrogen column density 
A^h ~ 10 22 cm -2 . The gas number density is n ~ 10 11 cm -3 , below the estimate by Ferland 
et al. (1996) and Baldwin et al. (1996) of n ~ 10 12 ' 5 cm -3 . We also show results for stellar 
atmospheres at d = 0.125 and 0.5 hg 5 pc, corresponding to half and twice the heating distance, and 
obtain similar results. Our results are insensitive to changing the "boundary" hydrogen number 
density by a factor of ~ 10. In Figure 6, we compare the f^ya that we derive with that which 
CLOUDY calculates; note that the escape probability approximation used by CLOUDY tends to 
overestimate the radiation pressure, at least in the cases examined here. 

As a check, we compare the excitation temperature derived from CLOUDY's escape 
probability approximation with that derived from our own calculation using equation (|L]). The 
results, shown in Figure 7, are quite consistent, implying that the equilibrium level populations 
are well described by the escape probability approximation. Hence, it is the scattering in the 
damping wings of the line profile, which does not significantly affect the level population but 
which greatly affects the line width, which is inaccurately calculated by CLOUDY in the escape 
probability approximation. In comparing the excitation temperature to the kinetic temperature, 
we see from Figure 7 that there is a region of column density ~ 10 22 cm -2 in which T exc exceeds 
Tkin by a factor ;> 2. Hence, the naive estimate of the radiation pressure using equation (||) may 
be as much factor of ~ 16 too small. 

Our objective has been to carry out an exploratory investigation of the importance of Lya 
pressure in the atmosphere of a star which is being irradiated by an AGN. This study has 
involved a number of simplifying approximations. In particular, we have treated the atmosphere 
as plane-parallel with normally-incident irradiation; this, of course, precludes examination of the 
important streaming effects which appear likely to ensue from the transverse pressure gradients 
which must arise from the non-uniform irradiation of the stellar surface. 

In view of the neglect of the transverse pressure gradients and resulting motions, we have 
not iterated to a fully self-consistent solution. Our hydrostatic equilibrium is obtained using 
CLOUDY's estimate for the radiation pressure rather than our more accurate result for Ph ya - Nor 
have we addressed the issue of stability or included radiation pressure from lines other than Lya 
and from the star's own luminosity. 

With these numerous simplifications, we find that in low gravity stars, the differential pressure 
due to Lya can rise to ~ 10% of the hydrostatic pressure. 
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5. Summary and Discussion 

We have developed a formalism for computing the Lya pressure in a plane-parallel stellar 
atmosphere. The treatment of resonance-line trapping given in §2 can be incorporated into already 
developed stellar atmosphere codes. For conditions believed representative of the BLR (ionizing 
flux of ~ 4 x 10 10 erg cm -2 s -1 ), stars with g <J 10 cm s -2 will develop an appreciable pressure due 
to trapped Lya photons in their atmospheres. 

We have neglected numerous complicating factors, all of which appear likely to enhance the 
rate of mass loss. First of all, the gravitational acceleration g oc r~ 2 , rather than the assumed 
g = constant. As a result, the outer layers of the atmosphere will have a larger scale height and 
must ultimately go over into a thermal wind, even if the lower layers are essentially hydrostatic 
and stable. With the outer layers of the atmosphere photoionized and heated to T4 ^> 1, a fluid 
element will have positive enthalpy at r ^ 600 T^ 1 (M/Mq)Rq. 

Furthermore, we have neglected the radiation pressure due to resonance lines other than 
Lya, as well as the radiation from the star itself. Ferland (1997) reports that including radiation 
pressure from all resonance lines rather than just a select few sometimes increases the total 
radiation by an order of magnitude (using CLOUDY's escape probability approximation, of 
course). Since no part of our derivation in §2 depends explicitly on the resonance line being Lya, 
it is straight forward to extend our calculation to include other resonance lines. 

Most important, however, is the fact that except along the "axis" — the line from the center 
of the star to the AGN — the atmosphere will be subject to substantial tangential stresses due to 
three effects related to anisotropic irradiation by the AGN: 

1. The momentum deposited by absorbed photons — at anywhere other than the axis, the 
momentum will have a tangential component; 

2. Transverse gradients in the Lya pressure — at any particular radial column density Ah, the 
Lya production rate will be greater along the axis, where X-rays are entering radially; and 

3. Transverse gradients in the gas temperature — the rate of X-ray heating will be largest 
along the axis, where the zenith angle is zero. 

These tangential stresses must result in significant tangential flows, with flow velocities that 
might approach escape speeds along the "terminator," perhaps resulting in a "cometary star" as 
envisioned by Edwards (1980) and Matthews (1983). 

It would obviously be of great interest to calculate axisymmetric fluid-dynamical models of 
red giant atmospheres subject to intense X-ray irradiation in order to investigate the effects of 
tangential stresses. It should be noted, of course, that we lack a quantitative theory for mass loss 
even from unperturbed evolved stars, so we are hardly poised to develop a definitive treatment 
including this new complication. Nevertheless, it does appear likely that trapped Lya may help 
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drive mass loss in low-g stars subject to intense X-ray irradiation. The extra Lya radiation 
pressure and associated tangential stresses might induce such evolving stars to "bloat" further and 
lose mass prematurely, so that the population of stars near the AGN core could be skewed from 
the average throughout the galaxy 

Finally we note that none of our calculations are dependent on the gravitating object being 
a star. Recently, Walker and Wardle (1998) have suggested the existence of a population of cool, 
self-gravitating clouds in the galactic halo to account for "Extreme Scattering Events" which 
occur in monitoring the flux of compact radio sources. They propose clouds of approximately a 
Jovian mass with a radius of about an AU, which would correspond to a gravitational acceleration 
g ~ 6 x 1CP 4 cm s~ 2 at the "surface." With a total column density of ~ 10 26 cm -2 , the majority 
of the AGN energy flux would be absorbed well before reaching (and disrupting) the center of the 
cloud. The proposed particle number density of ~ 10 12 cm -3 corresponds well to that inferred by 
Ferland et al. (1996) and Baldwin et al. (1996). Thus if these self-gravitating, sub-stellar clouds 
exist in host galaxies of AGN, they seem to be potential candidates for the origin of the BLR. It is 
unclear, though, how they could exist with a high metal content, as required by BLR observations, 
without rapidly cooling and collapsing (in much less than a Hubble time) prior to exposure to 
AGN irradiation. 

Although we conclude that X-ray induced Lya pressure would not be significant for main 
sequence stars near AGN, it could be important either for lower gravity stars which have evolved 
off the main sequence up the giant branch, or for substellar, self-gravitating clouds. It remains 
to be seen whether the resulting cometary structures can account for the required BLR covering 
factor and cloud population, and the AGN mass supply. 

We gratefully acknowledge Bohdan Paczyhski for useful and inspiring discussions, Gary 
Ferland for help and advice with CLOUDY, and R. H. Lupton for the availability of the SM 
plotting package. The research was supported in part by NSF grant AST-9619429. 
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Fig. 1. — The Lya photon injection rate per H nucleon as a function of total hydrogen column 
density A%, in the cases of 5 = 1, 10, and 100 cms -2 . The incident flux is given by equation ( |35| ) 
with 74/(47T(i 2 ) = 2.37 x 10 10 erg cm -2 s _1 , corresponding to a distance 0.25 /tg 5 pc from an AGN 
of luminosity L46 = 27.8 /ig 5 2 . The heating distance dh for this value of L46 is 0.26 ^ 5 pc. The 
atmosphere is assumed to have ran = 10 10 cm -3 at Ah = 0. The stellar radius -R*,ioo i n units of 
100i?Q is given for each value of g, assuming a stellar mass of Mq. 
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Fig. 2. — The Lya photon destruction probability per scattering eo as a function of total hydrogen 
column density Nn, calculated from equation (p9|), in the cases of g = 1, 10, and 100 cm s~ 2 . The 
incident flux is given by equation (|35| ) with A/(4ird 2 ) = 2.37 x 10 10 erg cm -2 s~ 1 . 



-17- 



10 3 




10 20 1Q 21 

N H (cm- 2 ; 



10 24 



Fig. 3. — The differential radiation pressure -P^ya (thick lines) due to trapped Lya photons is 
plotted as a function of total H column density iVn, for g = 1, 10, and 100 cms -2 . The incident 
flux is given by equation d35| ) with A/(Aird 2 ) = 2.37 x 10 10 erg cm" 2 s _1 . The thin lines are the 
pressures required for hydrostatic support P^yd at the same gravities [see equation (33)], where we 
have assumed a density n H = 10 10 cm at Nh = 0. 
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Fig. 6. — Comparison of the Lya pressure derived from our calculation with that calculated 
by CLOUDY in the escape probability approximation, for g = 10 cms -2 , and A/(4ird 2 ) = 
2.37 x 10 10 erg cm -2 s -1 as in Figures 1-3. 
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Fig. 7. — Comparison of the excitation temperature T cxc in Lya derived from our calculation 
with that calculated by CLOUDY in the escape probability approximation, for g = 10 cm s -2 , and 
A/(4-Kd 2 ) = 2.37 x 10 10 ergcm- 2 s -1 as in Figures 1-3 and 6. For reference, the kinetic temperature 
Tkm is also given. 



